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Estudio preliminar de la diversidad genética de Ilex guayusa ecuatoriana mediante marcadores inter 
simple sequence repeat (ISSR) 


Resumen 


Evaluamos el grado de diversidad genética de Ilex guayusa, una especie de importancia etnobotanica y comercial para las 
comunidades indigenas de la Amazonia ecuatoriana. Caracterizamos genéticamente a 157 individuos, provenientes de 
chacras ubicadas a lo largo de seis provincias, con nueve marcadores ISSR (Inter Simple Sequence Repeats). Se detecté en 
total 91 bandas polimorficas, pero el indice de heterocigosidad estimada (H, =0.19) de la poblacion reveld un nivel reducido 
de variabilidad genética para la especie. Los indices de distancias genéticas de Nei (0.013 < Ds < 0.086) revelaron un nivel 
reducido de divergencia genética entre individuos provenientes de provincias diferentes. De manera similar, el andlisis de 
varianza molecular (AMOVA) demostré que solamente el 18% de la variacidén observada para J. guayusa ocurre entre 
poblaciones. El bajo grado de diversidad genética encontrado para I.guayusa puede ser atribuido al hecho de que esta 
especie es cultivada exclusivamente mediante propagacién clonal; una actividad cultural que probablemente ha 
homogenizado el acervo genético de la especie a lo largo de su rango geografico de cultivo. El andlisis de coordenadas 
principales (PCoA) demostré que el germoplasma colectado puede estructurarse en tres grupos distintos caracterizados por 
un leve contraste genético en gradiente direccional, de norte a sur. La inclusi6bn de un mayor numero de muestras de 
provincias sub-representadas (ej. Sucumbios) y poblaciones silvestres, si existen, ayudaria a esclarecer el nivel de diversidad 


genética y estructura poblacional de la especie, y la historia de su cultivo en la regién. 


Palabras Clave. Aquifoliaceae, Guayusa, propagacién clonal, baja diversidad genética. 
Abstract 


We assessed the degree of genetic diversity of /lex guayusa, a species of ethnobotanic and commercial relevance for in- 
digenous communities of the Ecuadorian Amazon. We characterized 157 individuals, from small cultivation sites across six 
provinces, using nine Inter Simple Sequence Repeat (ISSR) markers. A total of 91 polymorphic bands were detected in, but 
estimated heterozygosity (H~ = 0.19) revealed a reduced level of genetic variability. Pairwise-Nei genetic distance indices 
(0.013 < Ds < 0.086) implied a reduced level of genetic divergence between individuals from different provinces. Ac- 
cordingly, partitioning of genetic diversity (AMOVA) indicated that only 18% of variation observed for . guayusa occurred 
between populations. The low degree of genetic diversity found for /. guayusa could be ascribed to the fact that the species 
is exclusively cultivated via clonal propagation; a cultural activity which has probably homogenized the species' genetic pool 
across its geographic range of cultivation. However, PCoA analysis resolved collected germplasm into three distinct groups 
displaying a subtle genetic contrast in a directional gradient, from north to south. The inclusion of a higher number of sam- 
ples from underrepresented provinces (i.e. Sucumbios), and wild populations, if they exist, would help resolve knowledge 


gaps regarding the genetic diversity and population structure of the species and its cultivation history in the region. 
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INTRODUCTION 


Ilex guayusa Loes. (Aquifoliaceae) is a tree species 
native to the Amazonian rainforests of northern South 
America [1]. In Ecuador, it is almost exclusively found 
in small cultivation sites (commonly known as chacras) 
along the western basin of the Amazon. It currently rep- 
resents the most harvested and widely utilized non-food 
plant species by indigenous groups, like the Kichwa, 
Achuar and Shuar [2,3]. These ethnic groups typically 
use the leaves of J. guayusa to prepare ceremonial bev- 
erages conferring good luck for hunting and fishing and 
for protection against snake bites [2-5]. Ilex guayusa 
leaves are also employed for ritualistic purging, and to 
treat ailments like gastritis, stress and infertility [5,6]. 
The leaves of J. guayusa are rich in caffeine, terpenes, 
phenols and other methylxanthines, and exhibit putative 
anti-inflammatory, microbicidal, anti-oxidant and anti- 
obesity properties [4,5,7]. 

Historically, dried 1. guayusa leaves were only sold in 
low volumes as a specialty product in local markets 
throughout the Ecuadorian Amazon. However, since 
2009 and principally owing to the efforts of RUNA?™, 
a social enterprise based in the USA and Ecuador, /. 
guayusa has gained commercial notoriety in international 
markets. It is commercialized as soluble teas and canned 
drinks with energizing and stimulant properties, mainly 
in the United States and Canada. New export markets 
have thus generated an increase in the number of I. gua- 
yusa producers in the Ecuadorian Amazon, from ap- 
proximately 300 in 2011, to over 1,600 in 2013 [8]. Cer- 
tainly, the growth in popularity of . guayusa can bring 
socio-economic change to the region, but further ex- 
pansion in its commercial production must be carefully 
evaluated as the species is native to a fragile and highly 
biodiverse ecosystem, vulnerable to the detrimental ef- 
fects of intensive monocropping [8]. 

Despite its rich ethnobotanic relevance and rising com- 
mercial prospects, the biology, cultivation history and 
genetics of Ilex guayusa are poorly understood. This 
research presents a first assessment of the genetic di- 
versity of I. guayusa in the Ecuadorian Amazon using 
Inter-Simple Sequence Repeat (ISSR) markers. This 
preliminary survey of the degree and geographic distri- 
bution of the genetic diversity of the species is neces- 
sary for establishment of diversity conservation strate- 
gies, design of plant improvement programs and gener- 
ation of sustainable-agroforestry systems that best uti- 
lize plant genetic resources in vulnerable ecosystems. 


MATERIALS AND METHODS 
Plant material and DNA isolation 


RUNA's technicians, in collaboration with students from 
USFQ, collected leaf samples from 157 Ilex guayusa 
individuals, from small cultivation sites listed in their 
database of J. guayusa leaf providers. The sites were 
distributed across the provinces of Sucumbios (1), Napo 
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(91), Orellana (7), Pastaza (35), Morona Santiago (11) 
and Zamora Chinchipe (12); all located in the Ecuado- 
rian Amazon (Fig 1). Collected samples were georef- 
erenced using a Garmin E-Trex Legend HCx GPS sys- 
tem (Garmin International Inc., USA). For every indi- 
vidual, young and healthy leaves were collected from 
the lower branches of trees and preserved at -20°C. Ge- 
nomic DNA was isolated from leaf samples using the 
CTAB procedure described by Kieleczawa [9]. Follow- 
ing extraction, DNA quality and concentration were as- 
sessed using a NanoDrop 1000 Spectrophotometer 
(Thermo Scientific, USA). DNA samples were diluted 
to a final concentration of 20 ng/yul. 


PCR amplification 


Nine anchored ISSR primers were used for the molec- 
ular characterization of Ilex guayusa (Table 1). These 
primers belong to the universal UBC ISSR primer set 
(University of British Columbia, Canada) and had been 
successfully employed by Zhou [10] for the genetic char- 
acterization of diverse Ilex species. Genomic DNA (40 
ng) was PCR amplified in a 20 jul reaction containing 50 
mM KCI, 20 mM TrisHCl (pH 8.4), 1 ug/l BSA, 2.0 
mM of MgClo, 0.16 mM Primer Mix, 0.5 mM dNTPs, 
and 2 units Taq Polymerase (Life Technologies, Carls- 
bad, California). PCR amplification conditions consisted 
of an initial denaturation at 94°C for 4 min; followed 
by 40 cycles of a 1 min denaturation at 94°C, anneal- 
ing of 50 sec at 52°C, an extension at 72°C for 1 min; 
and a final extension period of 7 min at 72°C. All reac- 
tions were performed in a T-Personal Series Thermocy- 
cler (Biometra, Goettingen, Germany). Amplified prod- 
ucts were separated and visualized via electrophoresis 
in 1.75% w/v agarose gels, run at 90 V for 2 hours, and 
stained with SYBR Safe DNA gel stain (Life Technolo- 
gies, Carlsbad, California). Gel band sizes were esti- 
mated using a | Kb DNA ladder (Life Technologies, 
Carlsbad, California) and band scoring was performed 
visually based on band-size differentiation. All bands 
were scored as presence/absence markers, and a binary 
raw-data matrix, where the presence of a band was en- 
coded as one (1) and the absence of a band as zero (0), 
was generated for all gel banding patterns. 


Data analysis 


For all identified loci expected heterozygosity (.) was 
determined using the GenAIEx 6.501 software package 
[11]. Accordingly, within- and between-population es- 
timates of molecular variance (AMOVA) and Nei ge- 
netic diversity indices were obtained using GenAlEx 
6.501; assuming individual provinces as independent 
genetic groups, and excluding Sucumbios because only 
one sample was available for this province. Principal 
coordinates analysis (PCoA) was performed with DAR- 
win 5.0.15 using the DICE genetic dissimilarity coeffi- 
cient [12]. 
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Figure 1: Location of the 157 I. guayusa individuals collected marked as pink circles. 


RESULTS 


Collectively, all evaluated ISSR markers yielded a to- 
tal 91 polymorphic bands across the complete sample 
set (157 individuals). The total number of polymorphic 
loci scored per primer ranged from 3 (UBC-860) to 16 
(UBC-857), with an average of 10.1 loci per primer. Of 
the 91 bands analyzed, 73 polymorphic loci occurred 
in low frequencies (<0.5%) across the sample-set, but 
none were exclusive to any given province or genotype 
group. The expected heterozygosity estimate for [lex 


guayusa (H, = 0.19) revealed a low degree of genetic 
diversity for the collected germplasm. Accordingly, pair- 
wise Nei genetic distance indices implied a reduced level 
of genetic variability for this species (Table 2). The 
greatest degree of genetic differentiation was detected 
between individuals from Napo and Zamora Chinchipe 
(D, = 0.086); the two most geographically distant pro- 
vinces. The lowest genetic distance was detected be- 
tween individuals from the neighboring provinces of Napo 
and Orellana (D, = 0.013). While a Mantel test did not 
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Table 1: Information of nine ISSR anchored primers used for the molecular characterization of 157 Ilex guayusa 


individuals from the Ecuadorian Amazon. 


Primer Primer Number of 
code sequence amplified bands 
UBC-815 = (CT)8G 11 
UBC-824 (TC)8G 7 
UBC-825 = (AC)8T 15 
UBC-826 (AC)8C 12 
UBC-827. = (AC)8G 15 
UBC-840 (GA)8YT 11 
UBC-846 (CA)8RT 11 
UBC-857  (AC)8YG 18 
UBC-860 (TG)8RA 6 
Total 106 


Number of Polymorphism 
polymorphic bands percentage 
10 90.91 
6 85.71 
13 86.67 
12 100 
14 93.33 
10 90.91 
7 63.64 
16 88.89 
3 50 
91 85.85 


Table 2: Nei genetic distance matrix comparing Ilex guayusa individuals from different provinces in the Ecuadorian 
Amazon. Sucumbios was excluded from the analysis as only one sample was available for this province. 


Napo Orellana’ Pastaza Morona Santiago 
0 
0.013 0 
0.04 0.032 0 
0.04 0.023 0.025 0 
0.086 0.075 0.044 0.027 


reveal a significant association between genetic and ge- 
ographic distances (r=0.060; p=0.037), geographically 
distant provinces generally displayed a higher degree of 
genetic divergence. 


Analysis of molecular variance (AMOVA) demonstrated 
that the highest percentage (82%) of genetic variation 
for Ilex guayusa was found within groups (individuals 
belonging to the same province), relative to among groups 
(18%). This factor indicates a low degree of popula- 
tion differentiation. Nevertheless, Principle Coordinates 


Analysis (PCoA) (Fig 2) could resolve collected germplasm 


into three robust genetic clusters as determined by its 
first two coordinates, each respectively accounting for 
23% and 11% of observed variation. Cluster A included 
individuals from the most-northern areas of the provinces 
of Napo and Orellana, as well as the only individual 
from the province of Sucumbios. Cluster B included in- 
dividuals from central and northern quadrants of Napo, 
and individuals from the province of Orellana that were 
geographically proximal to Napo's eastern provincial bor- 
der. Finally, Cluster C comprised all individuals from 
the southern quadrants of Napo, and individuals from 
the southern Amazonian provinces of Pastaza, Morona- 
Santiago and Zamora Chinchipe. 
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DISCUSSION 


The main objective of this study was to assess the level 
of genetic diversity of Ilex guayusa, a species of great 
ethnobotanic and commercial importance for indigenous 
communities of the Ecuadorian Amazon [2]. Overall, 
expected heterozygosity (H, = 0.19) revealed a reduced 
level of genetic variability for the species. For compar- 
ison, using five isoenzyme markers, Winge et al. [13] 
detected a comparatively higher degree of heterozygos- 
ity (H, = 0.50) for Ilex paraguariensis, a close relative 
of I. guayusa. The RAPD-based characterization of four 
natural populations of I. paraguariensis from southern 
Brazil revealed that the average (Nei) genetic distance 
between individuals within a population was 0.392 [14], 
whereas the maximum degree of genetic differentiation 
between J. guayusa individuals analyzed in this study 
was 0.086 (Table 2). At first hand, the low level of ge- 
netic variability observed for Ecuadorian J. guayusa ap- 
pears contradictory, especially when considering the ex- 
tensive number of polymorphic loci (91) detected across 
the evaluated sample-set. However, it is important to 
highlight that approximately 80% of all detected poly- 
morphisms occurred in low frequencies (<5%), and the 
vast majority of these did not exhibit correlative pat- 
terns of allelic fixation for specific genotypes or geno- 
type groups (i.e. individuals from specific provinces). 


The low degree of genetic diversity found for Ilex guayusa 
could most likely be ascribed to the fact that all col- 
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Figure 2: PCoA of the 157 Ilex guayusa individuals. Each number corresponds to an individual, from one of the 
six provinces of Ecuadorian Amazon, grouped in clusters A, B and C: Sucumbios (green), Napo (black), Orellana 
(orange), Pastaza (blue), Morona-Santiago (red) and Zamora-Chinchipe (purple). 


lected germplasm originates from small cultivated plots. 
In the Ecuadorian Amazon, indigenous groups cultivate 
I. guayusa exclusively via clonal propagation through 
stem cuttings, a cultural activity that has probably ho- 
mogenized the species’ genetic pool across its geographic 
range of cultivation [16-18]. Historically, . guayusa 
has existed in close relationship with Amazonian eth- 
nic and cultural processes, and the species' sampled ge- 
ographic distribution likely reflects the migration pat- 
terns of indigenous groups [15]. Moreover, no reports 
exist of wild /. guayusa populations, and it is therefore 
impractical to extrapolate on whether the species dis- 
plays a generalized low degree of genetic diversity. The 
results of this study also imply that in the Ecuadorian 
Amazon, cultivated 1. guayusa originates from a spe- 
cific and reduced gene pool. Partitioning of genetic di- 
versity (AMOVA) indicated that 82% of the variation 
observed for this species occurred within populations 
(provinces), and only 18% between populations. Natu- 
ral populations of I. paraguariensis have also shown a 
reduced degree of genetic differentiation, principally ac- 
counted to the species’ outcrossing nature [14]. In out- 
crossing perennial tree species, natural populations gen- 
erally exhibit a high degree of allelic exchange via gene 
flow; a factor which limits genetic divergence between 


populations [19]. Yet, because samples analyzed in this 
study come from cultivated plots established through 
clonal propagation, it seems probable that the high de- 
gree of genetic homogeneity observed for I. guayusa 
results from the widespread propagation of a selection 
of homogenous genotypes with desirable characteristics 
and established for cultivation prior to the anthropogenic- 
driven dispersion of the species across the Amazon basin 


[5]. 


Remarkably, while genetic differentiation between es- 
tablished groups was limited, a subtle genetic structure 
could be established for Ecuadorian Ilex guayusa. Prin- 
ciple Coordinates Analysis (PCoA) could resolve col- 
lected germplasm into three distinct groups. All clus- 
ters included a high proportion of individuals from the 
province of Napo, which was over-represented in this 
study. However, each cluster was distinctly and respec- 
tively characterized by the presence of individuals from 
the remaining northern, central and southern provinces 
of the Ecuadorian Amazon (Fig 2). For decades, Napo 
has represented the center of economic activity in the 
Ecuadorian Amazon, and around 50% of its population 
is composed of migrants from neighboring provinces. It 
is probable that some of these migrants brought specific 
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variants of I. guayusa for cultivation [20]. Incidentally, 
PCoA clustering tended to group individuals from Napo 
in a directional gradient, from north to south, possibly 
highlighting the geographically structured introduction 
of novel J. guayusa alleles from distinct, yet proximal 
geographic locations in neighboring provinces. Certainly, 
the inclusion of a higher proportion of individuals from 
underrepresented provinces (i.e. Sucumbios), as well as 
wild populations (if these should exist), would help us 
better resolve the geographic distribution of I. guayusa 
genetic diversity along the Amazon basin of Ecuador, 
and formulate stronger conclusions regarding its culti- 
vation history in the region. This information will be 
crucial for the establishment of adequate conservation 
strategies aimed at preserving and promoting the use of 
unique sources of genetic variation which can increase 
the productivity and environmental resilience of culti- 
vated orchards under sustainable agro-forestry systems. 
Expectantly, the latter will enable for a controlled and 
sustainable expansion of [. guayusa cultivation in the 
Ecuadorian Amazon which can meet rapidly growing 
commercial demands, and reduce the need for intensive 
mono-cropping systems and their associated detrimen- 
tal effects on vulnerable ecosystems. 
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